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Chaotic dynamics of a nonlinear oscillator is considered in the semiclassical approxi- 
mation. The Loschmidt echo is calculated for a time scale which is of the power law in 
semiclassical parameter. It is shown that an exponential decay of the Loschmidt echo is due 
to a Lyapunov exponent and it has a pure classical nature. 
PACS numbers: 05.45Mt, 05.60Gg, 03.65Sq. 
Classical chaotic dynamics can be characterized by a Lyapunov exponent A. The quantized procedure 
stops the classical spread of stretching and folding due to the uncertainty principle and, as a result, breaks 
the applicability of semiclassical approximation. The corresponding breaking time was found in [1] to be: 

rr, = (l/A)ln(/o/n), (1) 

where Iq is a characteristic action. It indicates a fast (exponential) growth of quantum corrections to the 
classical dynamics due to chaos. Recently this result gained a renewed interest also related to the fidelity of 
wave functions or the Loschmidt echo. In the field of quantum chaos a question on stability of trajectories 
can be re-addressed to stability of wave functions respect to a small variation of a control parameter [3] 
called the fidelity of the wave functions, which is a measure of a quantum reversibility. It is also named 
"Loschmidt echo" [4], and it is known that on the time scale of the order of Tn it decays exponentially ~ e~-^* 
[4]. Different regimes of decay of the Loschmidt echo have been observed for i.g. integrable [5,6], chaotic 
[5,7,8] and pseudo-integrable [9] systems (see also references therein). 

In the present paper, we show that the Loschmidt echo in the nonlinear kicked oscillator decays exponen- 
tially due to the Lyapunov exponent in a, so-called, Lyapunov regime [4,8]. An analytical expression for 
this behavior in the framework of the semiclassical expansion is obtained. 

Dynamics of a nonlinear oscillator with the Hamiltonian Hq — huja^a + fi^ ^{o^aY is integrable [10]. Here 
uj and ^ are linear frequency and nonlinearity parameters correspondingly, while, for the chosen notation, 
annihilation and creation operators have the commutation rule [a, a^] = 1. There is nontrivial semiclassical 
expansion that leads to an appearance of so-called Z?~forms. These forms are determined as derivatives over 
the initial conditions, say, (a, a*) as 

where A{B) = {ip{a)\A{B)\ip{a)) is the average value of the operator. Therefore, the D-forms determine 
the local instability of a dynamical system. In the presence of a perturbation dynamics becomes chaotic and 
the exponential growth of the D-forms due to the Lyapunov exponents leads to the logarithmic breaking 
time of Eq. (1), known also as the Ehrenfest time. That result could be transparently seen in the coherent 
state basis [1,2,10], and it is independent of the choice of the initial basis of wave functions [11]. Role of 
the D-forms as a sensitivity of dynamical variables to the initial conditions or to the perturbation can be 
also re-addressed to the wave functions. More detailed consideration on the D-forms for different dynamical 
realizations in the nonlinear oscillator is considered separately elsewhere [12]. 

Originally, a question on sensitivity of wave functions to a variance was asked in [3] to characterize quantum 
chaos by the fidelity of the wave function 

M{t) = KV-olexpli J {n + Sn)dt/h}exp{~i J ndt/h}\^l;o)\^ . (3) 

It characterizes an evolution of the initial wave function V'o governed by the two slightly different Hamilto- 
nians Ti and the variational Hamiltonian Ti, + SH. M{t) was also referred to as the "Loschmidt echo" [4], 
where dynamics of the initial wave function due to Ti. after time t is reversed back to the initial state with 
the variational Hamiltonian TL -\- STi. The dynamical decay of the overlap is characterized by the Lyapunov 
exponent on some semiclasical time scale Tsd that also stands to be determined. We will consider analyti- 
cally the overlap function M{t) for chaotic dynamics of the nonlinear oscillator in the presence of a periodic 
perturbation of the form of 5-kicks with a period T and an amplitude e. The Hamiltonian of the system is 
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H = nLoa''a + n^l2{a''af -Ueia'' + a) ^ d{t - nT) = Ho + V . 



(4) 



For the initial wave function we choose the coherent states. At the moment t = it determined as an 

cigcnfunction of the annihilation operator: a{t = 0)\a{t = 0)) = a\a) = a\a.). It enables one to study their 
dynamical evolution analytically in the semiclassical limit. Evolution of a generating function was in detail 
considered in [13]. We borrow some details of that analysis with corresponding adaptation to the present 
consideration. The evolution operator is 



U{t)=exp<.—i / dr [ujo'a + hfi{a' a)"^ — eg(T){a^ + a)] 



(5) 



where exp means T-ordered exponential. Under this symbol all exponents commute. Therefore, one can use 
the following Stratonovich-Hubbard transform [15] for the exponential 



\ rr f\ f i \2 /rj.] /■ TT dX{T) i r*drA^(T)/4K -if* 

exp —iTijjiT / drya'a) /T = / I I e Jo ^ " • e Jo 
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drA(r)a^a 



(6) 



where we use that k = hjit and t/T t is a, number of kicks represented in the continuous form. We 
take into account that the harmonic oscillator, acting on the coherent state, changes its phase only, and the 
perturbation V acts as a shift operator. Namely, e.g. for the period one writes 



i€{a^ +a) -i<i>xiT-)a^ a\ 



e"-^- , = exp[i(e/2)(ae-''^^(i) + a*e''^^(^)]|ae-''^^(i)"*" + ie) , 

where (/.a(I) = /□ dT[ujT + \{T)] is a phase on the period 1. Finely, we obtain the wave function in the form 
of the following functional integral [13] 



= U{t)\a) = / n (t^A(T)/V4^) exp 



drA^(r)/4K 



X exp 



le / dTg{T){al{T)+ax{T))l2 



where 



a{t) = e-'^^Wa(t) = g-'-^^W 



a + ie [ dTg{T)e''^^^^^ 
Jo 



(7) 



(8) 



Mt)= I dT[ujT + X{T)]. (9) 

Jo 

The wave function (7) is the evolution of the initial ket- vector \a) due to the Hamiltonian H of (4). To 
obtain an echo, we reverse dynamics at the moment t back to t = with a random {e.g. Gaussian) time- 
dependent process rj{t) to add it to the linear frequency u). This variation of the linear frequency affects 
efBciently the chaos control parameter K, since the last is K = AeijT\a{t)\'^ [1]. Therefore, evolution of the 
initial bra-vector (aj is due to the variational Hamiltonian [14] H + dH, with the random time-dependent 
frequency: 

w ^ = w -h r]{t)/T . 
The wave function of the variational motion is 

%{t) = [U{v,t)\a)]''. (10) 
For simplicity, we consider the Loschmidt echo averaged over the Gaussian distribution 

f dT{v{T)-fjf /4a 





P[r,{t)]^P{v{t),f],a) 



exp 
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with the non-zero first r] and the second a moments. It reads 

Ma{t)= I lld7i{T)P[ri{T)]M{rj,t) , 

where 

Mir^,t)^\{^,{m{m' 

is the Loschmidt echo for a one fixed realization of t]{t) with < t <t. Denoting by 



(11) 



(12) 



Px = -i dTg{T)ax{T) , 
Jo 

we obtain the scalar product of the wave function in (12) in the form 

MM ^M = mm{t)) = y n ^M^^exp [^jyr (A?(r) - A^(r)) 



X exp 



(13) 



In the semiclassical limit k <C 1 this integral is strongly simplified and its evaluation is analytically tractable. 
After linear change of variables 



Ai = 2/x + Kv 12, A2 = 2/z — Kv 12 — 77 , 



(14) 



the phases are 



^ = y drfwT -l-Ai] = y rfr[a;T + 2/i-|-K;z//2], 
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y dT\u:T + A2 + (5] = y dT\ujT - 2/x + kj//2]. 



The Jacobian of the transformation is 2k. Taking into account these expressions for the phases we obtain 
from (8) in the semiclassical approximation that 



"A2(*)"Ai(i) - e/3Ai - e/3A2 ~ / dTv(T)\a(T)Y 

Jo 



(15) 



[Q!Ai - c^As] ~ / dTv{T)a{T) . 
Jo 



(16) 



Here a(r), defined in (8), is taken for v{t) = 0. It gives the definition of the classical action as I{t) = K\a{t)\'^ 
for K ^ and |ci(t)|^ — > 00. To carry out integration over z/(t) we use the following auxiliary expression [13] 



exp 



-(^ 



K^/2)\ f dTiy{T)a{T)\A = — I d'^e-2|«lV«exp [-i^i?er f dTv{T)a 

Jo \ T^l^ J V ^ Jo 



Substituting (15), (16) and (17) into (13), we obtain for the scalar product of the wave functions 



dii{T)dv{T) 
2^r 



exp 



— / dT (4M(r) - 7?(r)) {kv{t) - r]{T)) 



exp 



IK / dTv{T){\a{T)+i\^-\e) 



(17) 



(18) 



The functional integrations over z/(t) is exact and gives the (5-function in /i. Hence the integration over /i(r) 
is also exact. After these integrations we obtain from (18) 



M = — d'ie 

KIT ' 



-2|?lV« 



exp 



- / dTri{T)Ici{a,a* ,t) 



(19) 



where Ici{a, Oi* ,t) = I{loT — 2|^p, a + ^, a* + r), and wo neglect a small term of the order of rj^ in the 
exponential. Let us expand the last exponential in (19) in the Taylor series in ^ and Therefore, we have 



(20) 



exp / dTv{T%i{a, a*,T)\ = T(loT - 2|CP, « + C, a* + T) 

= Sm,n,; n!m!i! a(f T)' •^('^^' ^i ^*) ' ^'"^*"(-2|^p)' . 

Substituting (20) in (19) and taking into account that 

— / d2$e-2|€lV-^P^*9 = v'(«/2)P+9vW'5p,« , 
we obtain an expression for the scalar product in the form of the expansion in the semiclassical parameter k 



M=y ^".^y ^ C r,,^L, ^i^T,a,a*) ■ (-2)'(ac/2) 

n.l 



(21) 



It should be stressed that the strongest contribution to the sum (21) for the same orders of n is due to the 
derivatives over the initial conditions, namely due to the D-iovm: 



Z)(r)^i3(/,/) = i^^^("'"'"* 



da 



dI{T, a, a* 
da* 



„2Ar 



Therefore, we obtain, approximately. 



Mr^{t) !v exp 



dTr]{T)Ici{T) exp - J drr,{T) J dT'r,{T')D{I,i{T), I,i{t')) 



(22) 



(23) 



Finely, using this approximation, we obtain that the Loschmidt echo is determined by the following Gaussian 
integral 



M,{t)^ / lldr]{T)P[v{T)]\M{v,t)\^ . 



(24) 



In the case of narrow packet cr <C 1 , that might be roughly approximated by a 5-function -P[??(t)] 
S{t]{t) — ff), we obtain 



M„{t) oc exp 



exp 



^ 2At 



(25) 



A possibility of such super-exponential drop off is also mentioned in [5] . For an arbitrary Gaussian distribu- 
tion, we obtain that the Loschmidt echo decays exponentially like exp [— Af]. For instance, in the opposite 
case, when cr 1, using the Fourier transform 



I 



di 



An 



e-CV4g-«A(t) 



where A[t) = dTe^'^r]{T), we can calculate the functional integral over r/(r) and then over ^. Finely, we 
have 



27r+ — (e 



2At 



1) 



-1/2 



OC 



-At 



(26) 
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These decays of M„{t) in (25) and (26) are determined by the classical Lyapunov exponent A. This result 
has pure classical nature and is independent of the semiclassical parameter k, and it survives in the classical 

case when k = 0, as well. 

In conclusion, we presented an analytical evaluation of the Loschmidt echo M{t) by means of semiclassical 
expansion. It leads to some restriction on time of the validity of semiclassical description which is definitely 
not coincided with the quantum breaking time Tn in Eq. (1). We shown that M(t) has sense on this 
semiclassical time scale and M(t) decays exponentially due to the ZJ-form according Eqs. (23) and (25). 

The validity of Eqs. (25) and (26) for this time scale can be obtained from the semiclassical expansion 
(15) and (16) in exponential (13) that is the semiclassical expansion for the linear oscillator exposed to the 
quasi-random field A. Therefore, the Jacobian of the Hamiltonian flow is 



J = det\d[a{t), a* {t)]/d [da, a*] 2k sin (pij,{t) j u{T)dT + j v{T)dT 



(27) 



It reflects that the Liouville theorem does not valid for the averages [10] Since v{t) is the quasi-random field 
with the complex Gaussian distribution defined at the Stratonovich-Hubbard transform in (6), we able to 
evaluate the integral in Eq. (27) substituting a mean value of v. The first moment equals zero: (j^(t)) = 0, 
that is why we use the second moment equaled (/^^(t)) = IQi/n. Substituting the square root from the 
modulus of {i^^{t)) in (27), we obtain for the Jacobian 

J w 1 - 8K^/^tsm(j)^ + IGnt'^ . 

Therefore, the validity of the semiclassical expansion carried out in the exponential (13) and, consequently, 
validity of the exponential decay of M{t) due to the Lyapunov exponent defined in Eq. (26) is determined 
by the following semiclassica time scale 

t < Tscl - 1/4^1/2 . 

In the semiclassical limit, when /t ^ 1 the following inequality is definitely true: Tsd ^ tr. This expression 
also might be a possible explanation of the observation in [8] of the exponential decay of M{t) in the Lyapunov 
regime for a kicked rotor on times much longer than the Ehrenfest time tr. 
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by the Minerva Center of Nonlinear Physics of Complex Systems. 



[1] G.P. Barman and G.M. Zaslavsky, Physica (Amsterdam) 91 A, 450 (1978). 
[2] G.M. Zaslavsky, Phys. Rep. 80, 157 (1981). 

[3] A. Peres, Phys. Rev. A 30, 1610 (1984). 

[4] R.A. Jalabcrt and H.M. Pastawsky, Phys. Rev. Lett. 86, 2490 (2001). 

[5] B. Eckhardt, J. Phys. A: Math. Gen. 36, 371 (2003). 

[6] G. Benenti, G. Casati, and G. Veble, preprint nlin/0304032. 

[7] P.G. Silvestrov, J. Tworzydlo, and C.W.J Beenakker, Phys. Rev. E 67, 025204 (2003). 
[8] J. Vanicek and E. Hcllcr, preprint quant-ph/0302192. 

[9] Ph. Jacquod, I. Adagideh, and C.W.J. Beenakker, Europhys. Lett. 61, 729 (2003). 
[10] G.P. Berman, A.M. lomin, and G.M. Zaslavsky, Physica D 4, 113 (1981). 

[11] G.P. Berman and A.I. lomin, Phys. Lett. A 95, 79 (1983); Theor. Math. Phys. 77, 1197 (1988) [Teor. Mat. Fiz. 
[12] A. lomin, unpublished. 

[13] V.V. Sokolov, Tcor. Mat. Fiz., 61, 128 (1984). 

[14] We use here the variational Hamiltonian with the variance SH = 5{t)a^a/T instead of the perturbed one to not 

confuse it with the time dependent perturbation V. 
[15] R.L. Stratonovich, Dokl. Akad. nauk SSSR, 115 1097 (1957) [Sov. Phys. Dokl. 2, 416 (1958)]; J. Hubbard, Phys. 

Rev. Lett. 3, 77 (1958). For more details see also B. Miihlschlegel in Functional integration and its applications, 

ed. A.M. Arthurs (Claredon Press, Oxford, 1975) p. 124. 



5 



